Cell Research

Neuromedin B (NMB) and gastrin-releasing peptide (GRP), two bombesin analogs, are endogenous itch-specific neuropeptides that induce histaminergic and nonhistaminergic itch, respectively. Their functions are mediated by two G protein-coupled bombesin receptors, NMBR and GRPR. Here we present cryo-electron microscopy(cryo-EM) structures of G protein coupled NMBR and GRPR bound to NMB and GRP, respectively. The structures reveal that both bombesin receptors contain an extended and deep pocket to adopt NMB and GRP, with the conserved C-terminal motif of GH(F/L)M from both peptides to contact the toggle switch residues for the receptor activation. Together with mutational and functional data, our structures reveal the mechanism of ligand selectivity and specific activation of the bombesin receptors. These findings also pave the way to facilitate rational design of therapies targeting bombesin receptors for the treatment of pruritus.

Since SARS-CoV-2 Omicron variant (B.1.1.529) was reported in November 2021, it has quickly spread to many countries and outcompeted the globally dominant Delta variant in several countries. The Omicron variant contains the largest number of mutations to date, with 32 mutations located at spike (S) glycoprotein, which raised great concern for its enhanced viral fitness and immune escape[1-4]. In this study, we reported the crystal structure of the receptor binding domain (RBD) of Omicron variant S glycoprotein bound to human ACE2 at a resolution of 2.6 [A]. Structural comparison, molecular dynamics simulation and binding free energy calculation collectively identified four key mutations (S477N, G496S, Q498R and N501Y) for the enhanced binding of ACE2 by the Omicron RBD compared to the WT RBD. Representative states of the WT and Omicron RBD-ACE2 systems were identified by Markov State Model, which provides a dynamic explanation for the enhanced binding of Omicron RBD. The effects of the mutations in the RBD for antibody recognition were analyzed, especially for the S371L/S373P/S375F substitutions significantly changing the local conformation of the residing loop to deactivate several class IV neutralizing antibodies.

The high mortality of severe 2019 novel coronavirus disease (COVID-19) cases is mainly caused by acute respiratory distress syndrome (ARDS), which is characterized by increased permeability of the alveolar epithelial barriers, pulmonary edema and consequently inflammatory tissue damage. Some but not all patients showed full functional recovery after the devastating lung damage, and so far there is little knowledge about the lung repair process1. Here by analyzing the bronchoalveolar lavage fluid (BALF) of COVID-19 patients through single cell RNA-sequencing (scRNA-Seq), we found that in severe (or critical) cases, there is remarkable expansion of TM4SF1+ and KRT5+ lung progenitor cells. The two distinct populations of progenitor cells could play crucial roles in alveolar cell regeneration and epithelial barrier re-establishment, respectively. In order to understand the function of KRT5+ progenitors in vivo, we transplanted a single KRT5+ cell-derived cell population into damaged mouse lung. Time-course single-cell transcriptomic analysis showed that the transplanted KRT5+ progenitors could long-term engrafted into host lung and differentiate into HOPX+ OCLN+ alveolar barrier cell which restored the epithelial barrier and efficiently prevented inflammatory cell infiltration. Similar barrier cells were also identified in some COVID-19 patients with massive leukocyte infiltration. Altogether this work uncovered the mechanism that how various lung progenitor cells work in concert to prevent and replenish alveoli loss post severe SARS-CoV-2 infection.

Human trunk development, including spine and spinal cord organogenesis, is a coordinated, orderly, and interdependent process with spatiotemporal tissue patterning. However, the underlying cellular and molecular mechanisms remain largely unclear due to the lack of an effective model that can simulate the early development of human body axis. Here, we reported the long-term patterning and dynamic morphogenesis of human trunk through the formation of spine-spinal cord organoids (SSCOs) self-organized from three-dimensional culture of human PSC-derived neuromesodermal progenitors (NMPs). The SSCOs resembled the morphogenetic features of spine and spinal cord along the anterior-posterior axis, and showed the chondro-osteogenic and neural trajectories consistent with developmental dynamics of spine and spinal cord in gestational embryo through single-cell RNA sequencing (scRNA-seq). In addition, we identified a new HMMR+ bipotent cell population with self-renewal ability and neural/mesodermal competence but distinct from NMPs, which may be involved in trunk development and represent an invaluable tool for disease modeling of spine- and spinal cord-related disorders. Graphic Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/549829v1_ufig1.gif" ALT="Figure 1"> View larger version (80K): org.highwire.dtl.DTLVardef@14bfd30org.highwire.dtl.DTLVardef@1db7fe8org.highwire.dtl.DTLVardef@1aa66a1org.highwire.dtl.DTLVardef@380932_HPS_FORMAT_FIGEXP M_FIG C_FIG

Bats are a major "viral reservoir" in nature and there is a great interest in not only the cell biology of their innate and adaptive immune systems, but also in the expression patterns of receptors used for cellular entry by viruses with potential cross-species transmission. To address this and other questions, we created a single-cell transcriptomic atlas of the Chinese horseshoe bat (Rhinolophus sinicus) which comprises 82,924 cells from 19 organs and tissues. This atlas provides a molecular characterization of numerous cell types from a variety of anatomical sites, and we used it to identify clusters of transcription features that define cell types across all of the surveyed organs. Analysis of viral entry receptor genes for known zoonotic viruses showed cell distribution patterns similar to that of humans, with higher expression levels in bat intestine epithelial cells. In terms of the immune system, CD8+ T cells are in high proportion with tissue-resident memory T cells, and long-lived effector memory nature killer (NK) T-like cells (KLRG1, GZMA and ITGA4 genes) are broadly distributed across the organs. Isolated lung primary bat pulmonary fibroblast (BPF) cells were used to evaluate innate immunity, and they showed a weak response to interferon {beta} and tumor necrosis factor- compared to their human counterparts, consistent with our transcriptional analysis. This compendium of transcriptome data provides a molecular foundation for understanding the cell identities, functions and cellular receptor characteristics for viral reservoirs and zoonotic transmission.

At the end of 2019, the SARS-CoV-2 induces an ongoing outbreak of pneumonia in China1, even more spread than SARS-CoV infection2. The entry of SARS-CoV into host cells mainly depends on the cell receptor (ACE2) recognition and spike protein cleavage-induced cell membrane fusion3,4. The spike protein of SARS-CoV-2 also binds to ACE2 with a similar affinity, whereas its spike protein cleavage remains unclear5,6. Here we show that an insertion sequence in the spike protein of SARS-CoV-2 enhances the cleavage efficiency, and besides pulmonary alveoli, intestinal and esophagus epithelium were also the target tissues of SARS-CoV-2. Compared with SARS-CoV, we found a SPRR insertion in the S1/S2 protease cleavage sites of SARS-CoV-2 spike protein increasing the cleavage efficiency by the protein sequence aligment and furin score calculation. Additionally, the insertion sequence facilitates the formation of an extended loop which was more suitable for protease recognition by the homology modeling and molicular docking. Furthermore, the single-cell transcriptomes identified that ACE2 and TMPRSSs are highly coexpressed in AT2 cells of lung, along with esophageal upper epithelial cells and absorptive enterocytes. Our results provide the bioinformatics evidence for the increased spike protein cleavage of SARS-CoV-2 and indicate its potential target cells.

IgG exists mainly as monomer, but recent studies suggest that IgG forms hexamer to mediate antibody functions. Although many structures of IgG monomer and its fragments are determined, there is no IgG hexamer structure at atomic level due to the weak Fc-Fc interactions. Here we engineered a hexameric IgG with IgM tailpiece fusion and determined the structure by cryo-EM. IgG-Fc hexamer forms hexagon symmetry with the six Fcs lie in a plane with Fc-Fc interaction in a mutual lock-and-key mode. This structure provides structural insights into Fc-Fc interaction of IgG and reveals molecular basis for its function.

Human Slack (hSlack) channel is abundantly expressed in the central nervous system and regulates numerous physiological processes. Here, we identified hSlack channel is a voltage dependent physical force induced desensitization potassium channel. In the cell-attached mode, hSlack exhibited large voltage dependent potassium currents without significant delayed rectification or C-type inactivation. After breaking the cell membrane resulting the inside-out mode, the potassium currents were largely inhibited and it transited to the strong voltage dependent delayed rectified state. Structural results suggested that the "S6 site" confers the voltage dependent delayed rectification of hSlack channel, like our recent finding in hEag2 channel (In revision). Strikingly, the physical force induced desensitization and gating kinetics transition are mediated by the gatekeeper lipids in the central pore. Our results reported a novel physical force sensing channel through the pore lipids mediating physical force induced desensitization and gating kinetics transition. It opens a new door to understand the molecular basis of mechanobiology.

Voltage gradient is a general physical cue that regulates diverse biological function through voltage-gated ion channels. How voltage sensing mediates ion flows remains unknown at the molecular level. Here, we report six conformations of the human Eag2 (hEag2) ranging from closed, pre-open, open, and pore dilation but non-conducting states captured by cryo-electron microscopy (cryo-EM). These multiple states illuminate dynamics of selectivity filter and ion permeating pathway with delayed rectifier property and Cole-Moore effect at the atomic level. Mechanistically, a short S4-S5 linker is coupled with the constrict sites to mediate voltage transducing in a non-domain-swapped configuration, resulting transitions for constrict sites of F464s and Q472s from gating to open state stabilizing for voltage energy transduction. Meanwhile, an additional ion occupied at positions S6 potassium ion confers the delayed rectifier property and Cole-Moore effects. These results provide novel insight into voltage transducing and potassium current across membrane, and shed light on the long-sought Cole-Moore effects.

Phylogenetically distant coronaviruses have evolved to use ACE2 as their common receptors, including NL63 and many Severe acute respiratory syndrome (SARS) coronavirus-related viruses. We recently reported two Middle East respiratory syndrome coronavirus (MERS-CoV) closely related bat merbecoviruses, NeoCoV and PDF-2180, use Angiotensin-converting enzyme 2 (ACE2) for entry. However, their host range and cross-species transmissibility remain unknown. Here, we characterized their species-specific receptor preference by testing ACE2 orthologs from 49 bats and 53 non-bat mammals. Both viruses exhibited broad receptor recognition spectra and are unable to use ACE2 orthologs from 24 species, mainly Yinpterochiropteran bats. Comparative analyses of bat ACE2 orthologs underscored four crucial host range determinants, all confirmed by subsequent functional assays in human and bat cells. Among them, residue 305, participating in a critical interaction, plays a crucial role in host tropism determination. NeoCoV-T510F, a mutation that enhances human ACE2 recognition, further expanded the potential host range via tighter interaction with an evolutionary conserved hydrophobic pocket. Our results elucidated the molecular basis for the species-specific ACE2 usage of MERS-related viruses across mammals and shed light on their zoonotic risks.

TACAN is not a mechanosensitive ion channel but significantly linked to the mechanical hyperalgesia. In this study, we show that the human TACAN is a homodimer with each monomer consisting of a body, a spring and a blade domains. The body domain contains six transmembrane helices that forms an independent channel. The spring domain adapts a loop-helix-loop configuration with the helix running within and parallel to the membrane. The blade domain is composed of two cytoplasmic helices. In addition, we found that all the helices of the body and the spring domains are specifically associated with membrane lipids. Particularly, a lipid core, residing within a cavity formed by the two body and spring domains, contacts with the helices from the body and spring domains and extends to reach two symmetrically arranged lipid clusters. These results extremely imply that the membrane lipids coordinate with the membrane-embedded protein to sense and transduce the mechanic signal.

Sterol O-acyltransferase 1 (SOAT1) is an endoplasmic reticulum (ER) resident, multi-transmembrane enzyme that belongs to the membrane-bound O-acyltransferase (MBOAT) family 1. It catalyzes the esterification of cholesterol to generate cholesteryl esters for cholesterol storage 2. SOAT1 is a target to treat several human diseases 3. However, its structure and mechanism remain elusive since its discovery. Here, we report the structure of human SOAT1 (hSOAT1) determined by cryo-EM. hSOAT1 is a tetramer consisted of a dimer of dimer. The structure of hSOAT1 dimer at 3.5 [A] resolution reveals that the small molecule inhibitor CI-976 binds inside the catalytic chamber and blocks the accessibility of the active site residues H460, N421 and W420. Our results pave the way for future mechanistic study and rational drug design of SOAT1 and other mammalian MBOAT family members.

The neuropeptide 26RFa, a member of the RF-amide peptide family, activates the pyroglutamylated RF-amide peptide receptor (QRFPR), a class A GPCR. The 26RFa/QRFPR system plays critical roles in energy homeostasis, making QRFPR an attractive drug target for treating obesity, diabetes, and eating disorders. However, the lack of structural information has hindered our understanding of the peptide recognition and regulatory mechanism of QRFPR, impeding drug design efforts. In this study, we determined the cryo-EM structure of the Gq-coupled QRFPR bound to 26RFa. The structure reveals a unique assembly mode of the receptor extracellular regions and the peptide N-terminus and elucidates the recognition mechanism of the C-terminal heptapeptide of 26RFa within the transmembrane binding pocket of QRFPR. The study also clarifies the similarities and distinctions in the binding pattern of the RF-amide moiety in five RF-amide peptides and the RY-amide segment in neuropeptide Y. These findings deepen our understanding of the RF-amide peptides recognition, aiding in the rational design of drugs targeting QRFPR and other RF-amide receptors.

The pandemic of COVID-19 caused by SARS-CoV-2 continues to spread around the world. Mutant strains of SARS-CoV-2 are constantly emerging. At present, Omicron variants have become mainstream. In this work, we carried out a systematic and comprehensive analysis of the reported spike protein antibodies, counting the antibodies epitopes and genotypes. We further comprehensively analyzed the impact of Omicron mutations on antibody epitopes and classified these antibodies according to their binding patterns. We found that the epitopes of one class of antibodies were significantly less affected by Omicron mutations than other classes. Binding and virus neutralization experiments show that such antibodies can effectively inhibit the immune escape of Omicron. Cryo-EM results show that this class of antibodies utilizes a conserved mechanism to neutralize SARS-CoV-2. Our results greatly help us deeply understand the impact of Omicron mutations. At the same time, it also provides guidance and insights for developing Omicron antibodies and vaccines.

Adenosine receptors, expressed across various tissues, play pivotal roles in physiological processes and are implicated in diverse diseases, including neurological disorders and inflammation, highlighting the therapeutic potential of receptor-selective agents. The Adenosine A3 receptor (A3R), the last identified adenosine receptor, is also activated by breakdown products of post-transcriptionally modified tRNA and exhibits dual roles in neuron, heart, and immune cells, and is often overexpressed in tumors, making it a target for anticancer therapy. Despite extensive studies on the other adenosine receptors, the structure and activation mechanism of A3R, especially by selective agonists like N6-methyladenosine (m6A) and namodenoson, remained elusive. Here, we identified N6-isopentenyl adenosine (i6A), a novel A3R-selective ligand, via comprehensive modified adenosine library screening. Cryo-EM analyses of A3R-Gi signaling complexes with two nonselective and three selective agonists revealed the structural basis for A3R activation. We further conducted structure-guided engineering of m6A-insensitive A3R, which would greatly facilitate future discoveries of the physiological functions of the selective activation of A3R by modified adenosines. Our results clarify the selective activation of adenosine receptors, providing the basis for future drug discovery.

Channelrhodopsins harvest the light and convert photons to the cellular ion flow. The ion selectivity and activation mechanism at the atomic level remains unknown. Here we describe cryo-EM structures for H. catenoides kalium channelrhodopsin (HcKCR1), its paralog, sodium selective channelrhodopsin (HcCCR), an open state of HcKCR1 (C110T), the voltage-dependent inwardly rectifier (D116N) and higher potassium selective channelrhodopsin (B1ChR2) from Bilabrum sp, illuminating the ion selectivity and activation mechanism. Briefly, the hourglass shaped lumen is occupied by the stepwise dehydrated potassium in both intracellular and extracellular side. The aromatic amino acids likely function as partial dehydrated potassium filter in the extracellular lumen, and intracellular dehydrated ion occupying layer chooses the right size of dehydrated ion, thus specifying ion selectivity and the higher dehydration capacity, the higher potassium selectivity. Furthermore, structural comparison of HcKCR1 and C110T suggested that the conformational changes of retinal triggers the extracellular side of TM6 extension as well as the retinal interaction residues motion, which then leads to ion flow. Our results not only uncovered the ion selectivity mechanism of potassium or sodium selective channelrhodopsins, but also elucidated their activation mechanism. It may provide a framework for designing next generation optogenetic tools.

Processing bodies (P-bodies) are the membraneless organelles that play critical roles in RNA storage and decay. Abnormalities in P-bodies contribute to diseases and developmental disorders. Huge efforts have been applied to the identification of the protein components in P-bodies, however, the dynamics of RNA components of P-bodies in human embryonic stem cells (hESCs) during maintenance and differentiation remain largely elusive. Here, we characterized the RNA profiles of P-bodies from HEK293T cells, hESCs, and hESC-derived mesodermal cells. The number of P-bodies decreases upon hESC differentiation towards mesodermal fate, accompanied with the decreased RNAs within P-bodies. By functional analysis of the P-body enriched and P-body depleted genes across different cell types, we discovered the cell type-specific enrichment of P-body-genes and the potential association with human diseases. We also captured the non-coding RNAs, including long intergenic non-coding RNAs and transposable elements in P-bodies in a cell type-specific manner. Furthermore, we verified the involvement P-bodies in regulating the differentiation of hESCs towards mesoderm by over-expression of LSM14A and knocking down of SPTAN1. In summary, we characterized the mRNAs and non-coding RNAs in P-bodies of hESCs and hESC-derived mesodermal cells, discovering the potential roles that P-bodies play for the precise differentiation of hESCs.

The spatial and temporal atlas of gene expression in the human embryo at early gestation is critical in understanding embryo development, organogenesis, and disease origins. We obtained the spatiotemporal transcriptome from 90 sagittal sections of 16 whole human embryos from 3 to 8 post-conception weeks by Stereo-seq with high resolution and ultra-large field, establishing the development trajectory/regulatory profiling of 49 organs. We uncovered the organ-specific regulons as potential lineage-determining factors and identified the new regulatory networks during heart and brain development. The atlas refines the key organs/cell types vulnerable to virus infection and genetic disorders, and, reveals the dynamics of allelic gene expression in specific organs at different stages. These results present the first comprehensive delineation of the spatiotemporal transcriptomic dynamics of human organogenesis. One Sentence SummaryThe spatiotemporal transcriptome atlas presents a comprehensive delineation of human embryogenesis after gastrulation.

The SARS-CoV-2 Omicron variant exhibits striking immune evasion and is spreading globally at an unprecedented speed. Understanding the underlying structural basis of the high transmissibility and greatly enhanced immune evasion of Omicron is of high importance. Here through cryo-EM analysis, we present both the closed and open states of the Omicron spike, which appear more compact than the counterparts of the G614 strain, potentially related to the Omicron substitution induced enhanced protomer-protomer and S1-S2 interactions. The closed state showing dominant population may indicate a conformational masking mechanism of immune evasion for Omicron spike. Moreover, we capture two states for the Omicron S/ACE2 complex with S binding one or two ACE2s, revealing that the substitutions on the Omicron RBM result in new salt bridges/H-bonds and more favorable electrostatic surface properties, together strengthened interaction with ACE2, in line with the higher ACE2 affinity of the Omicron relative to the G614 strain. Furthermore, we determine cryo-EM structures of the Omicron S/S3H3 Fab, an antibody able to cross-neutralize major variants of concern including Omicron, elucidating the structural basis for S3H3-mediated broad-spectrum neutralization. Our findings shed new lights on the high transmissibility and immune evasion of the Omicron variant and may also inform design of broadly effective vaccines against emerging variants.

PEGylated interferon- (PEGIFN) demonstrates promising therapeutic outcomes against chronic hepatitis B (CHB), whereas patient response to PEGIFN therapy remains unsatisfied. Shutdown of hepatitis B virus (HBV) antigens by RNA interference (RNAi) could enhance PEGIFN efficacy in CHB patients, whereas the underlying immunological mechanisms remain obscure. We performed studies by utilizing our newly established extracellular humanized IFNAR (IFNAR-hEC) mice. An in-house constructed small interfering RNAs (GalNac-siHBV) was administrated to mice either alone or in combination with PEGIFN. The phenotypic and functional characteristics of peripheral and organ-specific immune cells were assessed by flow cytometry, ELISpot, RNA sequencing (RNA-seq), and single-cell RNA-seq (scRNA-seq) analysis. Our results demonstrated that combined treatment with PEGIFN and RNAi exerted a synergistic and prolonged inhibition of HBsAg ([~]4log10 IU/mL, vs PBS) and induced a higher incidence of HBsAg seroconversion ([~]30%), comparing with either monotreatment. Mechanistically, combined therapy improved the functionality of global T and B cells, triggered increased anti-HBs producing B cells, and enhanced IFN{gamma}-producing T cells. scRNA-seq analysis revealed that the combined therapy reduced inhibitory B cell-B cell interaction, enhanced MHC-I signaling mediated T cell-T cell communication, and improved T cell-B cell crosstalk, thus improving the functionality of T and B cells. Enhanced MHC-II signaling networks across B cells and hepatocytes/Cd8+ T cells further promoted HBsAg seroconversion in the combined treatment groups. These results together provided scientific rationale and lessons for the combination of the two towards better therapeutic efficacy. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=133 SRC="FIGDIR/small/626539v1_ufig1.gif" ALT="Figure 1"> View larger version (35K): org.highwire.dtl.DTLVardef@84785aorg.highwire.dtl.DTLVardef@f984fborg.highwire.dtl.DTLVardef@189a966org.highwire.dtl.DTLVardef@d081f9_HPS_FORMAT_FIGEXP M_FIG C_FIG HighlightsO_LIShutting down HBsAg through RNA interference augmented the antiviral immune effects of PEGIFN in chronic HBV-carrier IFNAR-hEC mice. C_LIO_LICombined RNAi plus PEGIFN augmented the functionality of T cells, promoted B cell activation and class switch, but also exerted some suppressive effects on B cells. C_LIO_LIReduced inhibitory B cell-B cell interaction, enhanced MHC-I signaling between T cells and T cells, and improved T cell-B cell crosstalk, improved the functionality of T cells and B cells. C_LIO_LIEnhanced MHC-II signaling networks across B cells and hepatocytes/ Cd8+ T cells further promoted HBsAg seroconversion in RNAi plus PEGIFN combined treatment groups. C_LI